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THE   GENESIS   OF   LAME'S   EUATI 


By  Prof.  E.  H.  NEVILLE. 


THIS  note  was  provoked  by  the  assertion  in  the  chapter 
added  to  Modern  Analysis  in  the  third  edition,  that  "all 
the  analytical  proofs "  of  Lamp's  expression  for  Laplace's 
equation  in  terms  of  orthogonal  curvilinear  coordinates  "are 
extremely  long  and  cumbrous",  and  by  the  restriction  to  real 
variables  in  the  treatment  of  confocal  coordinates  by  \\  hittaker 
and  Watson. 

The  note  begins  therefore  with  an  analytical  transformation 
of  Laplace's  operator,  which  is  quite  short,  although  it  proceeds 
from  first  principles,  and  care  is  taken  to  explain  the  meaning 
of  the  steps  if  the  variables  are  complex. 

It  is  then  pointed  out  that  the  use  of  confocal  coordinates 
does  not  presuppose  the  variables  to  be  real,  and  the  familiar 
formulae  are  found  without  reference  to  geometry.  An  account 
of  the  origin  of  ellipsoidal  harmonics  and  of  the  relation  of  the 
different  species  to  Lamp's  equation  follows  naturally,  and  the 
note  ends  with  a  transformation  into  confocal  coordinates  of 
the  operator  introduced  by  Niven. 

There  is  no  novelty  in  the  results,  and  the  methods  are 
straightforward  ;  all  that  is  claimed  is  that  the  algebra  is  brief 
and  to  the  point. 


1.  If  real  variables  x,  y,  z,  the  coordinates  of  a  point  in 
space  referred  to  rectangular  axes,  are  functions  of  three  inde- 
pendent variables  *,,  ta,  *3,  the  direction  ratios  7",  jr\  y'3  of  a 
directed  tangent  to  the  curve  along  which  tr  varies  and  the 
other  two  variables  are  constant  have  the  values  dxjds,.,  di//dsr, 
dzjdsr,  where  sr  is  the  arc  of  the  directed  curve  measured 
from  some  tixed'point  upon  it.  If  dsrldtr  is  denoted  by  H(r), 
the  rate  of  change  of  any  function  0  along  the  curve  is  <f>rl  H(r), 
where  </>,.  denotes  90/9/y,  the  function  $  being  supposed  ex- 
pressed as  a  function  of  «„  tr  t3',  in  particular, 

1.1  •/*  =  */  Fr>,    7r'  =  */,/tf(,h    -r- 
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and  since 

1.2  2(7'Y=1, 

s 

it  follows  that 

1.3  #(,)•=£«;, 

where  with  Lame  we  write*  SF(a]  to  denote  F(x)  +F(y}  +F(z). 
The  transformation  from  x,y,  z  to  t^  tt,  ^3  is  triply  orthogonal  if 


1.4 
and  therefore  if 

1.5  /SV^u.  =  0, 

The  tangents  to  the  curves  of  reference  then  form  a  direct 
system  if  the  determinant  (7™),  whose  square  necessarily  has 
the  value  1,  itself  has  the  value  1,  not  —  1.  If  the  determinant 
is  equal  to  —  1,  an  interchange  of  t2  and  t3,  altering  only  the 
notation  of  the  transformation,  gives  a  determinant  with  the 
value  1,  and  the  same  effect  is  produced  by  changing  the  sign 
of  one  of  the  three  functions  ZT(]),  ZT(z),  /7(s)  ;  thus  without  any 
loss  of  generality  it  is  possible  to  assume  for  a  triply  orthogonal 
system 

1.6  (7")  =  1, 
that  is, 

1.7  9  (a,  y,  «)/9  K,  *„  0  =  H()}H<2)H(3). 

If  a  vector  has  components  £,  ij,  £  relative  to  the  axes  of 
coordinates  and  r,,  TS,  TS  relative  to  the  frame  composed  of  the 
tangents  to  the  curves  of  reference,  then  on  the  one  hand 


1.8  £ 

and  on  the  other  hand 

1.9  Tr  =  7rtf  +  7rIi?  +  7rtC,     ''=1,2,3; 

since  the   vector  is  arbitrary,  each   of  the  sets  of  formulae 
1.8,  1.9  is  equivalent  to  the  other. 

2.    When  a  transformation  is  made  from  a  set  of  three 
complex  variables  #,  y,  z  to  another  set  t},  t^  t3,  all  the  formula? 

*  Not  SF(x);  the  superiority  of  Lame's  notation  to  that  commonly  ascribed 
to  him  is  obvious  in  such  a  case  as  xFty  \  Sx-  4  y&$  /  oxoy  +  zo-(fr  ,'  dxdz,  which 
becomes  simply  <Sz«52<£  /  SxSu. 
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of  section  1  remain  true,  but  some  of  them  express  definitions. 
By  1.3  are  defined  the  squares  of  three  arc  functions;  in  this 
paper  we  suffer  no  loss  of  generality  by  supposing  that  not 
one  of  these  squares  is  identically  zero,  for  three  distinct 
directions  cannot  be  mutually  perpendicular  if  one  of  them  is 
nul  ;  we  can  therefore  suppose  also  that  the  range  of  the 
variable  is  so  restricted  as  to  render  the  arc  functions  them- 
selves single  valued.  The  arc  of  the  ?-th  "  curve  "  of  reference 
between  two  given  points  is  definable  as  fH(i.)dtr,  but  the 
curve  is  a  two-dimensional  region.  The  functions  defined  by 
1.1  are  called  direction  ratios,  and  they  satisfy  1.2,  which  is 
now  a  consequence  ot  1.3,  not  its  basis. 

The  transformation  from  x,  y,  z  to  £„  t2,  t3  is  called  triply 
orthogonal  if  1.5  is  true,  that  is,  if  1.4  holds,  and  the  trans- 
formation is  said  to  be  direct  if  in  addition  1.6  and  1.7  are 
satisfied  ;  there  is  no  loss  of  generality  in  supposing  the 
transformation  to  be  direct.  To  know  that  1.8  and  1.9  are 
still  equivalent  relations  between  any  two  sets  of  complex 
numbers  £,  ?;,  £  and  T,,  T2,  r3,  it  is  not  necessary  to  have 
a  theory  of  vectors  in  complex  space;  it  is  an  immediate 
algebraic  corollary  of  1.2  and  1.4  that  1.9  follows  from  1.8, 
and  1.6  is  sufficient  therefore  to  secure  that  1.8  follows  from 
1.9. 

From  the  set  of  formulas  1.1,  written  as 


combined  with  the  equalities 


we  have  a  set  of  important  relations  comprised  in  the  one 
typical  formula 

2.1  7."^,  +  7r'#(,),  =  7,"#(.)  +  7"#(.),, 

where  H(r)t,  H(t)r  denote  8#(,.)/9^  0#(,)/^,..     From  1.2,  with 
s,  t  for  r,  s, 

2.2 

t 

and  therefore  from  2.1 
2.3 


3.    For  any  function  $,  expressed  on  the  one  hand  as  a 
function  of  a',  ?/,  z,  and  on  the  other  as  a  function  of  /,,  t3,  £3, 


Prof.  Neville,  The  genesis  of  Lame's  equation.     341 
we  have 

31  J  JL      r-n          I        -1  I        JL      /•  1          O         9 

.1  m    —  oj  CC    ~r  ffl    ?y.-ra>.5,,        /   —  1,4,O, 

and  therefore,  from  1.1, 

3f\  ,      I  TT  rl  i        i        v2  »        i        r3  ,  1       o      O 

•  *  V-(r)^7    rj:'7    r      '7    T-'       J==l,-,  O. 

This  is  a  set  of  relations  of  the  form   1.9,  and  is  therefore 
equivalent   to 

3.3     <h  = 


If  ;//  denotes  0x,  a  double  application  of  3.3  gives 


that  is 

3.4  ^=  S  {  (  27  V/#<,> 


r        s 


and  by  addition  from  the  three  formulae  of  this  type  we  have 
at  once,  utilising  1.2,  1.4,  2.2  with  s  for  q,  and  2.3, 


87.       r 


since  identically 


where  n(r)'  denotes  H(l]H{s}B^IH(r},  the  formula  is  equiva- 
lent to 


r 

that  is,  to 
,6 

3-6 


the  familiar  result  discovered  by  Lame". 

4.  Let  a,  £,  c  be  three  fixed  finite  complex  numbers, 
subject  only  to  the  condition  that  no  two  of  them  are  equal. 
Then  if  x,  y,  z  are  complex  variables,  the  equation 

4.1  aj«/(a  +  e)+y'l  (b  +  0)  +  r/(c  +  6}  =  1, 

Z2 
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regarded  as  an  equation  in  a  variable  0,  lias  three  roots.  It 
x,  y,  z  are  unrestricted  tlie  roots  have  to  be  regarded  as  the 
three  branches  of  a  single  function,  but  if  the  joint  variation 
of  x,  y,  z  is  adequately  limited,  the  three  roots  are  separate 
single-valued  functions  of  x,  y,  z.  These  roots,  which  will  be 
denoted  by  #,,  #2,  #3,  form  the  set  of  confocal  variables  asso- 
ciated with  x,  y,  z  by  the  base  numbers  «,  b,  c. 

We  shall  write  a«,  be,  c»  for  «  +  0,  6  f  ^,  c  +  9  and  ap,  bp,  cp 
for  o  +  #p,  b  +  0p,  c  -f  8p  ;  obviously 


4.2 


In  using  Lame's  symbol  of  summation,  we  shall  write  SF  (c/,  n) 
for  F(a,  *}  +  F(b,  y}  +  F(c,  z}.  Thus  the  equation  4.1  can 
be  written 


For  any  value  of  Q  the  function 


will  be  denoted  by  9  (0). 

The  expression  of  x\  y'\  z*  as  explicit  functions  of  the 
confocal  variables  is  immediate,  for  (a  +  6}(b  +  6}  (c  +  6}  6  (#) 
is  a  cubic  polynomial  in  6  in  which  the  coefficient  of  bs  is  —  1, 
and  the  only  function  of  this  kind  with  0,,  02,  #3  for  roots  is 
-  &].  Hence  identically 

(^  -()}  (6  - 


and  decomposing  the  fraction  on  the  right  we  have 
u*      _         a,a,a3  J_    .   __M&__      J_ 

^J        t  .       —     /  /  \   /  \     *  /i      I       /  /  \   i  I  \  I      i 


and  therefore 
4.4    x2=- 


-b}(a  —  c}'lj    ~  (b  -  c}  (b  —a]  '  (c — a)  (c—b}' 

While  or,  ?/,  js"  are  seen  to  be  single-valued  functions  of 
0i?  6^  03,  the  variables  x,  y,  z  are  single-valued  only  if  re- 
strictions are  imposed  on  their  variation. 

Jt  is  assumed  thut  in  the  plane  associated  with  the  complex 
variable  B  there  are  cuts  from  the  three  points  -a,  —  b,  —  c 
to  infinity;  on  the  cut  plane  the  square  roots  of  the  functions 
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a  ,  b  ,  c  sire  severally  single-valued,  and  with  a  definite  choice 
among  these  square  roots  and  among  the  square  roots  ot  the 
constants  b  —  c,  c  —  «,  a  —  b,  the  functions  x,  y,  z  become 
uniform  within  their  restricted  domains. 

It  is  convenient  to  denote  the  product  of  selected  square 
roots  of  ap,  &p,  cp  by  Ap;  thus 

A  f  =  a  be  ,     «=  1,  2,  3, 
P        p  p  p'     * 

and  Ao  is  in  fact  uniform  through  the  range  considered. 

Differentiating   logarithmically    the   formula?*   in    4.4,    \ve 
have 


4.5  Xp  = 

and  therefore  with  the  notation  of  1.3 

4.6  Hv'=Snp'  =  iSU*ldp'. 

If  £,  y,  z  and  #,,  ^2,  ^3  are  regarded  as  constants, 


hence 

=  -  lim   Q  (W(6  -ffj 


from  4.3,  where  the  product  is  for  the  two  values  of  q  that  are 
different  from  p.  That  is  to  say,  7/(|)2,  #(2)2,  H(3*  are  given  in 
terms  of  the  coiifocal  variables  alone  by 

4.7 


From  4.2  we  have 

A.'/rfP  -  s«f  K  =  -  ft  -  ^)  ^'/^X, 

and  therefore  since  each  of  the  sums  on  the  left  has  the  value 
unity, 

4.8 


*  The  "  geometrical"  argument,  from  the  normal  to  the  surface  of  reference, 
is  more  illuminating,  but  I  am  deliberately  avoiding  any  use  of  complex  space. 
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whence  from  4.5 

4.9  Situ  =  Q 


tlie  transformation  to  con  focal  variables  is  triply  orthogonal, 
whatever  the  values  of  a,  b,  c. 

5.  Since  S(tl,  H(^,  H(z)  are  given  by  4.7  as  functions 
of  #,,  #2,  #3,  the  expression  of  Laplace's  operator  in  terms 
of  confocal  variables  is  deducible  from  3.6  ;  since  Ap  is  a 
function  of  6p  only,  the  formulae  in  4.7  imply  that  U(r)  ]  &rH(r) 
does  not  depend  on  6r,  whence 


5.1  —    —(^-  -i  )  = 


and  therefore 
5-2 


_a_ 

80, 


-  ^  V  V  =  (0,  ~  03)  A,  A, 


,  -    ,      a 
The  irrationality  is  only  apparent,  since 

53  A   A  (A    ^-A'^+ 

*"80A    >-80rr   A"   8^2 

or  explicitly 

5.4        Ar      -      A,-  -        =  aAc,-     0rr  + 


«,        Pr 

If  0  is  a  function  such  that 

where  ty{r)  depends   only  on    6(r},  the  equation   V" 
equivalent  to 

5.6        ft  -  0,)  *(1)  +  (8,  -  0.)  *(J)  +  (0,  -  0J  *w  =  0, 
that  is,  to 

*(,)     0t     1      =0, 
*M     6      1 

li'  2 

*w    e*    l 


=  0  is 
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and  therefore  there  are  multipliers  yl,  B  such  that 

5.7     *{1)  =  401  +  /?,     *u 
from  the  last  two  equations 


.showing  that  A  and  B  do  not  involve  #, :  similarly  A  and  B 
are  independent  of  #,  and  #3,  and  therefore  these  multipliers 
are  absolute  constants. 

6.  A  normal  solution  of  Laplace's  equation  in  confocal 
variables  is  a  solution  of  the  form  A  (#,)  M  (0,)  N  (0J,  the 
product  of  three  functions  each  of  a  single  variable.  Since 
x\  y\  z*  are  themselves  products  of  this  kind,  there  is  no 
greater  generality  in  dealing  with  a  function  a/yVAMN  than 
in  considering  AMN,  but  it  is  useful  to  have  the  constants 
?,  m,  n  at  our  disposal.  Writing 

6.1  <b  =  xlymz"i\.  M  N 

we  have 

$,.       K^       lxr       my,       nz . 

6.2  —  =  ~r  -\ \-  -    -  +  --'', 

*i  /i  /y*  it  & 

tp       j\r        a,          y          z 

where  K},  K2,  K3  denote  A(^,),  M  (^2),  N  (03),  and  primes 
indicate  differentiation  with  respect  to  the  variable  involved. 
Substituting  from  4.5, 

6.3  •"-      "" 
and  therefore 


ar        r      c,. 


since  from  6.3 


5.4  and  6.4  imply 

/    80  \ 

(  Ar  ^J  ) 

\       d^r/ 

(  I       m       n\  *        [I        m 

+  i  —  +  r  +  ~   ~  a  (  —  i  +  r^  + 

\«r      ^r      cr/          \«r       V 


a  ,  /      »»    w 

^-  +  —   +  T-  +    -      ^r 

Kr  \ar      br      cr)  Kr 

I       m       n  *  I        m        n 


1        1        IN  (KJ          !  I       m 


r      cr> 
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that  is 

"    -    " 


-  ' 

do1,./         Kr     [_  \  ar          b 

(fl      m      n\fl+l     m+1     n+l\        (I 

+  ii  J-+r+     r  •+-r:-+-    -)-2(- 


ar       ,.     c,\ar  r          cr  /        \a        b,.      cr 


Tt  follows  that  the  function  xlijmzn  A  (#,)  M  (#J  N  (03)  is  a 
solution  of  Laplace's  equation  if  and  only  if  the  functions 
A  (#),  M  (#),  N  (^)  are  three  solutions  of  the  differential  equation 


6.6      -j— ^  -f-  ( -       ^  +  7 —      •  +  -       -  ] 


C+0J  d0 


m  n    ]     f  I  I  1        m  4  1 

•     1      i     /i   T  7\  I        i  n     i      "7 77    T 


[a  +  0      b  +  0      c  +  B\      a+&       b  v  V 


where -4,  -S  are  arbitrary  constants;  with  a  modification  of 
these  constants,  the  equation  takes  the  form 


6  7 h    — — -    H — ^  + 

•  7/J2        'I  i      /J      '        L      ,      /J       ' 

aa         \a  +  u       b  +  t7       c  -t-  c//  ac/ 
^(7-1)      ni  (»»-!)  .,  tiQi-1)  A'6  +  B' 

i   4,  1       .    /i  2  "i       /;.    .    /i\2  /  _   .    /j\  a          /  __   .    «i\  /L   i    zii  /  _   i.  *j\  f 


To  multiply  each  of  the  functions  -ff"r(#r)  ^J  flrxbr'icr\ 
where  X,  yu,,  v  are  constants  independent  of  r,  is  lo  multiply 
the  product  AMN  by  a  constant  multiple  of  x'^y^z".  Hence 
if  K(0]  is  a  solution  of  6.6,  (a  +  B)x  (b  4  BY  (c+  0}v  K(6)  is  a 
solution  of  the  equation  derived  from  6.6  by  replacing  /,  m,  n 
by  I—  2X,,  m  —  2/z,  n  —  2v.  In  particular,  the  solutions  of  6.6 
are  the  solutions  of  the  equation 


d\K    LfJ_      J^      J.,, 

d6***\at0  +  b+0*  d±0)  dd  ~±(a+V)(b+B)(c±V) 

divided  by  (a  +  0/''2  (b  +  0/"/2  (c  +  ftf2.  This  equation  G.8, 
which  is  Lame's  equation,  corresponds  to  zero  values  of  ^,  w,  ?/, 
and  is  therefore  the  equation  satisfied  by  each  of  the  three 
factors  of  any  harmonic  function  that  is  normal  in  the  confocal 
variables. 
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7.  If  k  has  one  of  the  values  -a,  —  Z>,  —  c,  the  product 
(0,  —  A-)  (fft  —  k)  (0,  —  /i)  is  a  constant  multiple  of  #',//,  or  23 
as  the  ease  may  be;  for  any  other  value  of  k  this  product  is 
seen  by  4.3  to  have  the  value  (a  +  k)(b-\  k}(c  +  k]  9  (&),  that 
is,  to  be  a  multiple  of  Su*j(d-\  k)  —  I  by  a  factor  independent 
of  the  variables  a;,  y,  z. 

A  solution  0  of  La  pi  ace' a  equation  which  has  the  form 
given  by 

/    a*  ?/*  z1  \fc> 

7.1      d)  = 


where  6\  /,  ?»,  H,  Ict,  q^  k,f  q^  ...  are  constants  and  the  number 
of  factors  is  finite,  is  called  an  ellipsoidal  harmonic.     Since 


is  a  constant  multiple  of 

nn(*r,* 

^    r 

that  is,  ofZ(0l)Z(ei)Z(03),  where 
7.2 


the  function  ^(^)  defined  by  7.2  must  be  a  solution  of  the 
equation  6.6  in  order  that  the  function  0  given  by  7.1  may 
satisfy  Laplace's  equation.  Since  no  solution  of  6.6  can  have 
an  accessible  pole  or  branch  point  except  at  one  of  the  critical 
points  —a,  —b,  —  c,  and  since  the  constant  kp  cannot  have 
either  of  these  critical  values,  the  index  g  must  be  a  positive 
integer,  and  therefore  the  function  Z(ff]  must  be  a  polynomial 
in  6,  of  which  neither  a  +  6,  b  +  Q,  nor  c  +  6  is  a  factor. 
Conversely,  any  polynomial  which  satisfies  6.6,  resolved  into 
factors,  provides  an  ellipsoidal  harmonic  that  is  determinate 
but  for  a  constant  factor. 

It  is  not  to  be  expected  that  for  arbitrary  values  of  the 
constants  involved  the  equation  6.6  should  possess  an  integral 
that  is  a  mere  polynomial  in  6.  Let  Pt  (z)  denote  genericallv 
a  power  series  p^  +  p^  -+  pyz!t+...  in  which  p0  is  not  zero,  and 
Pni  (z)  a  power  series  in  which  neither  p0  nor  pt  is  zero.  Then 
if  in  Lame's  equation  we  substitute  («  +  0)XP0  (a  +  0}  for  K, 
the  lowest  degree  in  (a  +  6}  that  occurs  is  \  —  2,  and  the 
coefficient  of  (a  +  0)x"  is  ^0{\(X-  1)  +  i\},  that  is,  p0M\-i). 
Hence  one  integral  of  the  equation  is  expansible  near  —a  as 
P0(a  +  6}  and  one  integral  as  (a  +  6)*PQ  (a  +  6),  and  any  integral 
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which  is  not  a  multiple  of  one  of  these  two  is  expansible  as 
P9i  («  +  #)*.  It  follows  that  the  equation  6.6  has  an  integral 
A0  whose  form  near  —  a  is  (a  +  6]  '2P0  (a  +  6),  and  an  integral 
At  whose  form  near  —a  is  (a  +  6}~V~1)Ppll(a-t  6),  and  that 
every  integral  which  is  not  a  multiple  of  either  Aa  or  Ar  is  of 

the  form  (a  +  6)~ '  JP01  (a  +  #)*  near  —  a.  Hence  unless  each  of 
the  numbers  Z,  wz,  n  has  one  of  tlie  two  values  0,  1,  every 
solution  of  6.6  has  either  a  zero,  a  pole,  or  a  branch-point  at 
one  of  the  three  points  —  a,  —  5,  —  o;  in  particular,  a  poly- 
nomial solution  without  one  of  the  factors  a  +  $,  b  +  6,  c+6  is 
then  impossible.  Thus  the  values  of  I,  wi,  n  are  restricted,  not 
by  any  analogy  with  questions  in  spherical  harmonics,  but  by 
the  nature  of  the  problem  proposed  in  the  definition  of  an 
ellipsoidal  harmonic,  and  ellipsoidal  harmonics  fall  naturally 
into  four  species"  since  of  the  values  of  the  three  indices 
I,  in,  n,  all  may  be  zero,  two  zero  and  one  unity,  one  zero  and 
two  unity,  or  all  unity. 

When  one  of  the  eight  possible  sets  of  values  has  been 
chosen  for  I,  wi,  n,  the  values  of  A  and  B  have  to  be  considered. 
If  the  degree  of  the  polynomial  Z(&]  is  /?,  the  coefficient  of 
the  highest  power  of  d  when  Z(Q]  is  substituted  for  K  in 
6.6  is 

h(h-l)  +  h(l  +  m  +  n  +  f ) 

+  ±{(l  +  m  +  n}  (I -f  m  +  n  +  3)  -  2  (I  +  m  +  «)  -  A], 
and  for  this  to  vanish 

A  =  ^hl  +  2h  (2l+2in  +  2n  +  l)  +  (l+m  f  n]  (I  +  in  +  n  +  l), 
that  is, 

7.3  A=q(q  +  \), 
where 

7.4  q  =  2h  +  l  +  m  +  n, 

that  is,  where  q  is  the  degree  of  0  in  the  variables  x,  y,  z. 

Thus  we  are  led  to  regard  as  the  most  important  case  of 
Lame's  equation  the  case  in  which  A  is  of  the  form  q  (q  +  1), 
where  q  is  a  positive  integer,  and  it  is  of  this  case  that  complete 
solutions  have  been  given,  for  an  arbitrary  value  of  B,  by 
Hermite  and  Halphen. 
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But  an  arbitrary  value  of  B  does  not  permit  a  polynomial 
solution  of  6.6.     If  in  6.6,  multiplied  thoughout  by 


h 

we  substitute  K(6)  = 

0 

where  the  coefficients^,  /j,  ...,  fh  are  at  our  disposal,  the 
coefficient  of  6htl  vanishes  on  account  of  7.3  and  7.4,  the 
coefficients  of  0A,  6h~l  are  of  the  forms 


the  coefficient  of  #'  for  every  value  of  s  from  h  —  3  to  1  is  of 
the  form 


and  the  term  independent  of  6  is  of  the  form 


where         &0,  A,,  ...,  5A,  cc,  c,,  ...,  cft_,,  rf,,  <7,,  ...,  dh_9 

do  not  involve  B  and  are  polynomials  in  a,  5,  c,  with  coefficients 
depending  on  /,  ?»,  ?*,  h.  If  we  take^  to  be  unity  and  equate 
to  zero  the  coefficients  of  6h,  6h~l,  ...,  6'\  9,  we  obtain  values 
of  the  coefficients  /A_p  /A_3,  ...v/j,  /„  in  succession  as  polynomials 
in  a,  J,  c,  ^?,  and  the  coefficient  of  B3  in  /A_f  is 


Equating  finally  to  zero  the  term  independent  of  0,  we  have 
an  equation  of  degree  h+1  in  B,  in  which  the  coefficient  of 
Bh+l  is  a  constant  different  from  zero  and  the  coefficients  of  the 
other  powers  of  B  are  polynomials  in  a,  £,  c.  There  is  a 
polynomial  solution  of  6.6  of  degree  h  if  and  only  if  B  satisfies 
this  equation  ;  to  each  root  of  the  equation  in  B  corresponds 
one  and  only  one  polynomial,  and  different  values  of  B  lead 
necessarily  to  different  values  of^_,  and  therefore  to  different 
polynomials  in  6  and  to  different  ellipsoidal  harmonics.  Lame, 
Liouville,  and  others  have  shewn  that  if  a,  /?,  c  are  real  and 
distinct  the  roots  of  the  equation  in  B  are  also  real  and 
distinct,  from  which  it  follows  that  the  number  of  distinct 
harmonics  to  which  the  equation  leads  is  exactly  A  +  l;  hence 
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in  this  case  the  number  of  distinct  ellipsoidal  harmonics  of 
degree  q  is 

if  q  is  even  and  is 


if  ^  is  odd,  that  is,  is  always  2^  +  1,  and  these  are  real  functions 
of  a?,  ?/,  x.  For  complex  values  of  a,  3,  c,  it  is  possible,  but 
exceptional,  for  the  equation  in  B  to  have  equal  roots,  and 
therefore  for  the  number  of  ellipsoidal  harmonics  of  degree 
q  to  fall  below  2q  +  l. 

8.  Niven  has  shewn  that  ellipsoidal  harmonics  are  derivable 
from  harmonics  homogeneous  in  x,  y,  z  by  operators  con- 
structed from  the  operator  a  (d'/dx2)  +  b  (3*/9/)  +  c  (3"/9*')) 
that  is,  Sdcfldu3,  which  is  denoted  by  D\  It  is  interesting 
to  express  this  operator  in  terms  of  the  confocal  variables, 
although  there  is  no  obvious  use  for  the  transformation,  since 
it  is  essentially  to  functions  of  x,  ?/,  z  that  the  operators  have 
been  applied.  To  derive  an  expression  for  D~§  from  3.4,  all 
that  is  necessary  is  a  determination  of  the  two  sums 


for  all  values  of  r  and  s,  where  dl,  d'\  d3  are  used  as  alternative 
symbols  for  a,  b,  c. 

Writing  the  identity 


in  the  form  SJu'jdd  =1, 

that  is,  2 

t 

we  have  from  1.2 

8.1  S^(7^ 

and  from  1.4 

8.2  2 

By  differentiation  of  8.1 

8.3  2  <*W'  =  -  (H(p)p  +  2H(p)*)  1 
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Also  from  2.1 

HM  2  d  '7  V  =  HU  2  ^  V  V  +  #(.),  2  d'  (7'T  -  //(,,,  S  d  V  V  ', 

t  «  «  t 

that  is, 


1  (7 


-* 


rom  8.1,  and  therefore 
8.4      v^'          =  - 


if  s  is  different  from  ?•,  the  variable  6r  enters  into  H(j'  only  in 
the  factor  !/(#,.—  0s),  hence 


=1/2  (0r  -  tfj  —  ^(,,. 
and  8.4  can  be  written 
8.5 


We  have  now  from  3.4 


* 


, 
' 


(r>    0i-   I 

~J 


Jt'  \vc  write 
8.6 


the  first  group  of  terms  here  is  E'ty  and  the  second  is  E^. 
Because  £rH(r)  is  a  homogeneous  function  of  0{,  6^  0,  of  the 
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first  degree, 

vfl'-8(A,.ff(,.)) 

r-~~80T      '  0)' 

and  because  A  .  is  a  function  of  9r  only, 

v  fl  8A'-  _  o  3A,.  e 
,   '80/   S0,  ' 
hence 

8.7  AS 


and  so  finally 
8.3       ^  = 

that  is, 
8.9    i  (9,  - 


From  3.3  and  1.1, 


for  any  triply  orthogonal  transformation  ;  if  the  variables  are 
confocal,  4.5  implies 


Hence  for  confocal  variables 


and  ^  is  a  homogeneous  function  of  degree  n  in  the  rectangular 
coordinates  x,  y,  z  if  and  only  if  E$  =  n<f>. 
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